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Using a simple phenomenological parametrization of the reaction amplitude we 
calculated a-particle spectrum in the reaction p+^^B— > a +^ Be* 3a at the 
^ , resonance proton energy 675 KeV. The parametrization includes Breit-Wigner factor 

^ ' with an energy dependent width for intermediate ^Be* state and the Coulomb and 

the centrifugal factors in a-particle emission vertexes. The shape of the spectrum 
\ consists of a well defined peak corresponding to emission of the primary a and a flat 



shoulder going down to very low energy. We found that below 1.5 MeV there are 
17.5% of a's and below 1 MeV there are 11% of them. 
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. P. , I. INTRODUCTION 

Successful realization of aneutronic fusion idea would greatly reduce problems associated 
with neutron radiation such as ionizing damage, neutron activation, and requirements for 
biological shielding, remote handling, and safety issues. One of the reactions discussed in 
;his respect is the reaction p+^^B— > a Be* 3a. There are at least two proposals 
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11, 121] how to use this reaction for energy production. The shape of a-particles spectrum, 
especially its low energy tail is of crucial importance for realisation of these projects. Existing 
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measurements j^, 4| produce the spectra for fixed angles of a-particles that does not say 
much about integrated spectrum. In this paper we propose a simple theoretical model for 
the reaction amplitude allowing to calculate with some assumptions the shape of a-particle 
spectrum. 



II. STRUCTURE OF THE REACTION AMPLITUDE 

The reaction p+^^B^ 3a has a well pronounced resonance at the proton energy 675 KeV. 
The resonance corresponds to excited state of ^^C nucleus with the energy 16.57 MeV, and 
with quantum numbers 2~ and isospin T = 1. The dominant decay mode is a-particle plus 
excited state of ^Be with the energy E* = 3.06 MeV and quantum numbers 2"*" and T = 0. 
The cross section for decay into the ground state of ^Be is about two orders of magnitude 
lower |5l]. The decay proceeds with isospin violation, therefore the width of the resonance is 
small enough. With the Q- value Qi = 6.14 MeV the total width is T = 0.2 MeV. It leads to 
high cross section at the resonance peak ~ 1.2 b. The decaying state has negative parity, 
therefore the primary emitted a must be in the state with an odd angular momentum. In 
our case it can be L = 1, and L = 3. 

The width of the excited ^Be state is 1.513 MeV. It corresponds to lifetime 4.35 x 10^^^ 
s. During this time the primary a-particle is still within strong interaction range. It means 
that the decay is in fact a 3-body type. Immediate consequence of it is that the energy 
of any a-particle can vary from to 2Q/3, where Q is total Q-value for the reaction, 
Q = Mc + 16.57 - 3M„ = 9.295 MeV. 

Keeping in mind these features of the reaction the amplitude can be presented in the 
following general form: 

Ti = V2{P2 + Pi/2) G(-pi) l^i(pi), (1) 

where ^i(pi) is the emission vertex of the primary a-particle with the momentum pi, G is 
the propagator of the intermediate system a-particle plus ^Be*, and V2(p2 + Pi/2) is the 
decay vertex of ^Be* depending on relative momentum of two a-particles. The diagram 
corresponding to this process is shown in Fig.l. 

For unpolarized initial state an angular distribution of a-particle is isotropic. For this 
reason we omit angular dependence in the vertexes ^i(p) and V2(p) and will discuss only 
their energy dependence. There are two factors determining momentum dependence of the 
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FIG. 1: Diagram of the reaction. Dashed hnes correspond to emitted a-particles. 

vertexes at low momentum. First, it is a Coulomb barrier that suppresses emission of low 
energy particles. Second, it is a centrifugal barrier. With these factors the decay vertex of 
^Be* can be presented as 

2 

V2{p) = 92 J- exp (-27rQ;/^;)r(l + lia/v), (2) 

V + P 

where ^; is a relative velocity of a-particles in cm. frame, and r(a;) is the Euler Gamma- 
function. The value of the cut-off parameter K2 is ^2 ~ where R is the nuclear radius. 
The shape of the a-particle spectrum appeared insensitive to exact value of K2- With this 

vertex we can define the energy dependent width by equation 

Here a is the fine structure constant, and p — ^/m^, where is the mass of a-particle, 
and e is the cm. energy. Remaining parameter 92 is determined by the width of the excited 
state of ^Be 

r2(e*) = 1.513, (4) 
where e* = 3.06 MeV. The vertex Vi{p) has similar structure 

Vi{p) = gi ^ ^ exp {-^'Ka/v)V{l + i%a/v). (5) 

Here the Coulomb barrier is larger, ZiX Z2 — 8. As for centrifugal barrier, strictly speaking 
there are two partial waves with L — 1 and L = 3 contributing to this vertex. However, since 
we are interested in low energy tail of the a's spectrum, we retained only L — 1. In addition, 
the total width of the excited carbon state is small enough, Fi = 0.2 MeV, therefore, the 

momentum dependence is smooth in this energy range and does not influence the shape of 
the spectrum. As in the previous vertex the exact value of the cut-off parameter ki is not 
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important. The shape of the spectrum is insensitive to it. The parameter gi is defined by 
normahzation of the spectrum to one. 

A propagator G of the intermediate a +^ Be* state is of the form 

G = (6) 
E-H 

where E is the energy of the initial Carbon state E = Mc + 16.57 MeV, and H is the 
Hamiltonian of a +^ Be* system. Omitting a Be* interaction and retaining only the 
width of Be state we obtain 

"^^"^'^ ~ Mc + 16.57 - 3.06 - Mse - tp, - Ep, + tr2{e)/2 ' 

Here ep and Ep are the laboratory kinetic energies of a-particle and ^Be. Introducing Q- value 
for this transition we obtain, 

^^-P'^ = n 1 , ^ r r W9 ' 

Qi - 1.5ep, + tT2{e2)/2 

where we used the relation Ep = O.bep. The center of mass energy e in r2(e) is defined by 
total Q-value minus kinetic energy of a and ^Be, e = Q — 1.5e(pi). 

In the final state we have three identical Bose-particles. Therefore, the amplitude eq.([I]) 
should be properly symmetrized. The full amphtude is as follows 

y _ ^2(P2 + Pl/2)\4(pi) ^ l^2(Pl+P2/2)l^l(p2) ^ ^2(P2/2 - Pl/2) Vi{-pi - P2) 



Qi - 1.5ep, + zr2(e2)/2 Q, - 1.5ep, + zr2(ei)/2 Q, - 1.5ep,+p, + zr2(e2i)/2 

(9) 



III. SPECTRUM OF THE a-PARTICLES 

With the amplitude eq.Q we can calculate the spectrum of a-particles. 

dW = 27r6{Q - ep, - ep, - ^v^^^.W?'^^^ ■ (10) 

As usual, the line over the amplitude squared means averaging over initial spin projections 
and summation over angular momentum projections of a-particles. After that the spectrum 
becomes isotropic. Strictly speaking, the interference terms in the amplitude squared may 
depend on the angle between pi and p2. However, the scalar product (pi ■ P2) is fixed by 
delta- function in eq. (|TO|) . The delta-function in eq.( |TOi) means that we neglected the small 
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width of the initial excited Carbon state. Integration over angles removes delta-function, 
we obtain 

maPidpip2dp2 T7^mldeide2 

ITK — . (11) 



dW = |T|2 



47]-3 47r3 

Integrating eq. (fTT|) over €2 we obtain the final equation for a-particle spectrum dW/dE 



dW Z"'™" 



(12) 



where e™„ = {^Q/2-?,E/A - ^/E/2f, and e^ax = {^/Q/2-3E/4: + ^fE j^f . The nor- 
malization factor A defined by the condition 

io dE 



-dE = 1. 



The calculated spectrum is shown in Fig. 2. The spectrum has a well pronounced peak 
corresponding to emission of primary a-particle and a broad shoulder going down to very 
low energy. This low energy tail contains 11% of a-particles with the energy less than 1 
MeV, and 17,5% of a-particles with the energy less than 1.5 MeV. 

dW/dE 




E (MeV) 



FIG. 2: Normalized to one the spectrum of a-particles 



It is worth to see how important are the Coulomb and the centrifugal factors included 
in the vertexes in eq.([2]). To do this we calculated the spectrum using constant F2 = 1.513 
MeV and constant Vi and V2. In Fig. 3 we show the spectrum obtained in this way (dashed 
line) and the full spectrum from Fig. 2. The main difference is in the peak position which is 
shifted downward for ~ 0.5 MeV and an increase of number of a-particles with the energies 
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FIG. 3: Comparison of two spectra obtained for energy dependent (full line) and energy indepen- 
dent (dashed line) r2 and the vertexes 

between 2.5 and 3.8 MeV. Note that in this region all three a-particles have the energies 
close to each other. In this region the effect of interference of three terms in eq. ([9]) is strong. 
The low energy tail is almost insensitive to it. 

In summary, using a simple phenomenological model for the reaction amplitude the shape 
of the a-particle spectrum has been calculated. Basing on the obtained shape we estimated 
the number of low energy a's. We found 17.5% of a's below 1.5 MeV and 11% below 1 MeV. 
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